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Abstract 

The universal bifurcation property of the Henon map in parameter space is studied 
with symbolic dynamics. The universal-L region is defined to characterize the bifur- 
cation universality. It is found that the universal-L region for relative small L is not 
restricted to very small h values. These results show that it is also a universal phe- 
nomenon that universal sequences with short period can be found in many nonlinear 
dissipative systems. 



PACS number: 05.45. +b 



I. INTRODUCTION 

One of the standard ways of investigating the dynamics of physical systems is by exploiting 
the universal (system-independent) property^"^ of them. The best understood transition 
sequence is the period-doubling cascade, which has been observed in a variety of physical 
systems. Beyond the accumulation point for the period-doubling sequence there is chaos. 
Two decades ago Metropolis, Stein and Stein showed^ that there is an ordered sequence of 
distinct periodic windows, each of which occurs for some range of control parameter, within 
the chaotic region for unimodal maps, = /(/i,x„). They have called this sequence 

the U-sequence since the ordering of the windows is system independent. Remarkably, this 
universality is also observed in systems with many degrees of freedom both experimentally^'^'^ 
and theoretically^"^ although the phase portraits of these two- or high- dimensional system 
still exhibit very complex behaviour which is clearly not one-dimensional or close to one- 
dimensional. It has been found that the periodic windows interspersed in chaotic region for 
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these systems are ordered in a systematic way as those of one-dimensional (ID) maps. The 
most striking and detailed observation is obtained in the Lorenz equations 

X = 10{y-x), 

y — rx — xz — y, (1) 
z — xy — 8z/3. 

On the parameter r axis with 45 < r < 400, all of the found 68 periodic windows of the 
Lorenz equations can fit into those of a ID antisymmetrical map with only one exception^. 
Experimentally, even though the Belousov-Zhabotinskii reaction involves more than thirty 
chemical species, it exhibits rather complex bifurcation behaviour that is modeled well by 
ID maps^. 

Despite these numerical and experimental observations, the underlying mechanism for 
the universal property is not fully understood. The motivation of this paper is to present 
an approach towards interpreting all these experimental and numerical observations and 
exploring their limitations. We will take the Henon map ^'^ 

yn+l = bXn, 

as an example. The bifurcation structure of the Henon map in the two-dimensional pa- 
rameter (a, 6) space has been extensively discussed^^'^^. In this paper, we will use symbolic 
dynamics^'^^~^° of ID mappings and 2D mappings to illustrate the universal topological 
property of the Henon map at selected parameters by considering the unstable periodic or- 
bits embedded in its chaotic attractor. Two topological quantities 5 and L are defined to 
characterize this universal topological property. Then we discuss the universal bifurcation 
property of the Henon map in 2D parameter space (a, h) by defining universal-L regions in 
which the Henon map exhibits ID bifurcation behavior to period L. It is remarkable to 
find that the universal-L region for relative small L is not restricted to very small h values. 
Wc will also present two examples of ordinary differential equations (ODE's), the Rosslcr 
equations^'^ and the forced Brusselator'', to demonstrate the validity and robustness of our 
approach. These results show that it is also a universal phenomenon that universal sequences 
with short period can be found in many experiments or numerical calculations on nonlinear 
dissipative systems. 

The paper is organized as follows. In Section 11, we review the basic property of ID 
unimodal maps. The universal bifurcation property and its limitations of the Henon map in 
the 2D parameter (a, h) space is studied in Section III. To demonstrate the validity of the 
method presented in Section HI, the universal bifurcation property of Rosslcr equations and 
the forced Brusselator in a definite parameter axis is investigated in Section IV. Finally, in 
Section V we give our conclusion. 



II. UNIVERSAL SEQUENCES IN ID UNIMODAL MAPS 

By using the symbolic dynamics of ID mappings. Metropolis, Stein and Stein (MSS) had 
already shown that the dynamics of unimodal ID maps of the interval [-1, 1] is embodies 
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in the U-sequence of periodic windows Fig. 1 shows a typical case. The extremum 
is denoted by a letter C. Each periodic window of the map can be labelled by a symbolic 
sequence of O's and I's that mark the location (to the left or right of C ) of the successive 
iterates of the initial point C. For example, the only windows with period 5 are lOl^C, lO^lC, 
and lO^C ( lOl^C represents the periodic window (101^C)°^ hereafter.) Indeed, we can define 
an ordering^^'^^ for these symbolic sequences referring to the natural order in the interval 
[-1, 1]. These ordering rules are consistent with the ordering of a real number a defined for 
a sequence S(x) with an initial point x as foUowing^^ 

oo 

a{^{x)) = Y^ixa-\ (3) 

with [16] 

Since the symbolic sequence K=S(C), called also the kneading sequence, acquires a maximal 
a in this metric representation, a symbolic sequence S(x) corresponds to a real trajectory if 
and only if it satisfies 

Q;((7'"(S(a;))) < q;(K), m = 0, 1,2,---, (4) 

where a denotes the shift operator. With this admissibility condition, we can generate all the 
admissible periodic orbits for a given kneading sequence K. The kneading sequence changes 
as the controlling parameter alters. Since kneading sequences correspond to orbits coming 
from C, they should also satisfy the above condition. Thus we obtain the admissibility 
condition for K themselves: a symbolic sequence K can be a kneading sequence if and only 
if it satisfies 

q;((7'"(K)) < q;(K), m = 0,l,2,---. (5) 

When K is a periodic string, K corresponds to a periodic window. From Eq. 5, we can 
generate all the possible periodic windows. It can be checked that there are only three 
period 5 windows as those listed above. 

With the ordering rules in equation (3) , all periodic windows can be ordered to yield the 
U-sequence. In the logistic map, this U-sequence is consistent with the increasing order 
which is hsted in Table 1 up to period 7. 



III. UNIVERSAL SEQUENCES IN 2D HENON MAPS 

The Henon map (2) has be extensively studied by using symbolic dynamics^^""*^^. The set of 
all "primary" tangencies between stable and unstable manifolds determines a binary gener- 
ating partition which divides the attractor into two parts marked by letters and 1. Any 
trajectory is encoded by a bi-infinite string S(a;) = • • • • • • s_iSo • S1S2 • • • Sn " " " ) where s„ 

denotes a letter for the nth image, a letter for the mth preimage, each is either or 1, 
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Table 1: Symbolic sequences for periodic windows of the Henon map along two different 
parameter axes and that for the forced Brusselator equations. The axis I (long dashed in 
Fig. 4) and axis II (dash-dotted) are two axes in and out of the universal-7 region of the 
Henon map, on which a complete and incomplete U-sequence is found respectively. 



No. 


Period 


Word 


a range 


Brusselator* 








b 


=0 


Axis I 






Axis II 


w 








C 





-0.749 


1 


032565- 


1 


164538 


none 






2 


IC 





75-1.249 


1 


164539- 


1 


332034 


none 


0.45-0.544 




2x2 


lOlC 


1 


25-1.367 


1 


332035- 


1 


409865 


none 


0.545-0.5777 


2 


6 


lOlllC 


1 


4747-1.47973 


1 


499840- 


1 


495505 


none 


0.58249-0.58251 


3 


7 


lOllllC 


1 


57472-1.57541 


1 


585860- 


1 


585225 


1 


521515-1.522540 




4 


5 


lOllC 


1 


6244-1.62843 


1 


635395- 


1 


631640 


1 


588650-1.594165 


0.5845-0.5848 


5 


7 


lOllOlC 


1 


67396-1.6744 


1 


679020- 


1 


678610 


1 


651665-1.652290 




6 


3 


IOC 


1 


75-1.76853 


1 


763210- 


1 


743840 


1 


776860-1.792900 


0.5947-0.654 




2x3 


lOOlOC 


1 


76854-1.77722 


1 


772255- 


1 


763215 


1 


792920-1.800525 


0.6545-0.7025 


7 


7 


lOOlOlC 


1 


8323-1.83239 


1 


835474- 


1 


835567 


1 


811364-1.811581 




8 


5 


lOOlC 


1 


86059-1.86136 


1 


863083- 


1 


863841 


1 


843810-1.844705 


0.7068-0.7115 


9 


7 


lOOlllC 


1 


8848-1.88483 


1 


886869- 


1 


886911 


1 


870425-1.870470 




01 


6 


lOOllC 


1 


90726-1.90736 


1 


909005- 


1 


909120 


1 


894790-1.894915 


0.718-0.7185 


00 


7 


lOOllOC 


1 


92715-1.92716 


1 


925261- 


1 


925289 


1 


939390-1.939415 




12 


4 


lOOC 


1 


94056-1.94153 


1 


938829- 


1 


939827 


1 


952015-1.952945 


0.7345-0.792 


13 


7 


lOOOlOC 


1 


95371-1.95371 


1 


952122- 


1 


952138 


1 


964545-1.964555 




14 


6 


lOOOlC 


1 


96677-1.9668 


1 


968064- 


1 


968099 


1 


956760-1.956795 


none 


15 


7 


lOOOllC 


1 


977179-1.977184 


1 


978379- 


1 


978385 


1 


967795-1.967795 




16 


5 


lOOOC 


1 


98541-1.985468 


1 


984101- 


1 


984160 


none 


0.8259-0.8675 


17 


7 


lOOOOlC 


1 


991814-1.991818 


1 


992925- 


1 


992928 


1.982968-1.982969 




18 


6 


lOOOOC 


1 


996375-1.996379 


1 


995150- 


1 


995153 


none 


0.9015-0.923 


19 


7 


lOOOOOC 


1 


999096 


1 


997890- 


1 


997891 


none 





the solid dot indicates the "present" position. In order to extend the grammar for unimodal 
maps to this map, a "backward" variable is defined as^^ 

oo 

P{^{x)) = Y.i,a-^'+'\ (6) 

i=l 

with 

i., = I for I ^i=° Z ? (mod 2) for b < 0, 

For this 2D map, each primary tangency C is associated with a bi-infinitc kneading 
sequence K (with the first backward letter Sq undetermined which may be or 1) and two 
symmetrical points (q;(K), /3_ (K)) and {o.{K), /?+ (K) = 1 - /3_ (K)) in the symbolic plane 
corresponding to sq = and 1 respectively Analogously to those in unimodal maps, for all 
admissible points (a, /3) with (3 e (K), /3+ (K)], a should be less than a{K) and thus the 
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pruning front''^''' is obtained by cutting out rectangles {a,P\a > a{K),P e [P_{K), P+{K)]} 
for all points on the partition. The union of these rectangles gives fundamentally forbidden 
zone. Consequently, the grammar for a word admissible or forbidden in this map can be 
expressed as: A bi-infinite word is admissible if and only if all its shifts never fall into 
the fundamentally forbidden zone^'''^^. It is clear that there are infinitely many kneading 
sequences (corresponding to infinitely many primary tangencies) in a 2D map to determine 
the admissibility condition for a word, while there is only one kneading sequence in a ID 
map. 

Universality in the Hcnon map. Fig. 3 shows a typical symbolic plane, (a, 6) = (1.4, 
0.16). The corresponding attractor is shown in Fig. 2 which has a rather complicated 
structure. Its fractal dimension is 1.16± 0.03. Numerically 203 kneading sequences are 
found as shown in Fig. 2. It is found that the minimal and maximal of all the forward 
parts of these kneading sequences start with Kmm=101111010101 and Km(i^=101111011111 
respectively, corresponding to a minimal and maximal a-values a„jj„ = 0.837560 and cxmax 
=0.838466 of all these kneading sequences. We define two quantities 5 and L as 

5 = Oijnax - Oimin = 0.000906, (7) 

L = -[/o^25] = 10, (8) 

where [/o5f25] denotes the integer part of log25. It is clear that 5 =0 and L — >■ +00 in the ID 
hmit (6 =0). For (a, h) = (1.4,0.16), an unstable periodic orbit with length n < L = 10 can 
not tell the difference between these kneading sequences. Indeed, no symbolic string with 
length n < L lies in the interval between ^^min and Kmax- Thus for the unstable periodic 
orbits with length n < L, the grammar is completely determined by a symbolic string K/, 
which is 1011110101 or 1011110111, the first 10 letters of Kmin or K^ax, that is, a word S{x) 
corresponds to an unstable periodic orbit of the Henon map for (a, 6) = (1.4,0.16) if and 
only if it satisfies 

aia'^iSix))) <a{Kf), m = 0,l,2,---. (9) 

This is just the grammar for unimodal maps with a kneading sequence Kf. Consequently 
the unstable periodic orbits of the Henon map for (a, 6) = (1.4, 0.16) can be generated as 
that of unimodal maps with a kneading sequence Kf (see Eq. 4). The only exception is the 
unstable periodic orbit K^^ which can not be determined by Eq. 9. We here noted that the 
Henon map is divergent for a — 1.4 and b > 0.315. 

As b decreases, L increases. It had already shown^^ that L = 32 for (a, 6) = (1.4, 0.05). 
In the 2D phase space, even for (a, b) = (1.4, 0.05) the attractor has a clear hook indicating 
that the map is two-dimensional. We emphasize that though the attractors reveal very 
complicated structure in 2D phase space, the topologies for these attractors may be very 
close to those in ID maps that the unstable periodic orbits can be generated with only one 
kneading sequence to some degree. 

Now we consider the universal bifurcation property of the Henon map in parameter space. 
Fig. 4(a,b) show the isoperiodic lines^^'^^ for all the nine period 7 windows. Numerically we 
find that L < 7 for all the parameters a and b in the region between the two heavy solod 
lines shown in the figure^^. We call this region the universal — 7 region hereafter. Thus 
all the periodic orbits with length<7 of the Henon map in this Universal-7 region can be 
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determined with only one kneading sequences as those of ID unimodal maps. Consequently, 
in this region there is a perfect MSS-sequence up to period<7 along any axis provided that 
the axis is never tangent with any isoperiodic lines. These axes are in a sense the same as 
the axis of 6 = (corresponding to the Logistic map). We present an example of these axes 
in Fig. 4 (line 1, long dashed). The periodic windows on this axis are listed in Table 1. It 
is clear that they share the universal feature as that of ID unimodal maps up to period 7. 
We can also obtain universal — M region numerically for M = 5, 6, 8, 9, • • • in which there 
are MSS-sequences for period< M along any axes provided that they are never tangent with 
any isoperiodic lines for period< M. 

In Fig. 4 we also show the borders for the Henon map exhibit an attracting set with 
initial points (xo,?/o) very close to original point (0,0). Comparing to this borders we can 
say that the universal-7 region is not restricted to very small b values. Thus it is rather likely 
to get a MSS-sequence up to a relative short period (say, period 7) in the full 2D parameter 
plane of the Henon map. 

Incomplete U — sequence in the Henon map. In fact, even on a axis out of the uni- 
versal region, the Henon map can exhibit approximately ID behaviour if the axis is never 
tangent with any isoperiodic lines. In tabic 1 we also show the periodic windows on the 
axis represented by dash-dotted line (11) in Fig. 4. It is clear that all of these words in- 
crease monotonically as a increases except the word lOOOlC and the period windows lOlllC, 
lOOOC, lOOOOC and lOOOOOC are missing. 



IV. APPLICATIONS TO ODE's 

The above idea can be extended to many other two- or higher-dimensional systems. Here 
we only take the Rossler's equations^^ 

X — —y — X, 

y^x + ay, (10) 
z — b + z{x — c), 

and the forced Brussclator^ 

i; = A — (B + l)x — x^y + acos{u!t), 
y = Bx — x^y, 

as examples. 

The 2D attractor of the Rossler's equations is usually taken from a section of the 3D 
flow on the half-plane |/ = 0, x < 0^^. It has already shown that the unstable periodic 
orbits of the attractor can be generated with only one kneading sequence up to period 12 
for parameters c = 2, = 4 and a = 0.408 (corresponding to L > 12). We find similar 
results (L > 9) for c = 2, d = 4 and 0.125 < a < 0.415. Table 2 shows the periodic windows 
up to period 9 in descending a order along with their periods, words and locations on the 
parameter axis. They are exactly consistent with part of the U-sequence from word C up to 
lOOlOllC. 
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Table 2: The symbolic sequences for the periodic windows of the Rossler's equations 



No. 


Period 


Word 


a range 








C 


0.125-0.335 




2 


IC 


0.336-0.375 




2x2 


lOlC 


0.376-0.3834 




2^x2 


lOlllOlC 


0.3836-0.3852 


2 


6 


lOlllC 


0.390668-0.39097 


3 


8 


lOlllllC 


0.393624-0.393638 


4 


9 


lOllllllC 


0.395446-0.395449 


5 


7 


lOllllC 


0.396638-0.396676 


6 


9 


lOllllOlC 


0.398034-0.398041 


7 


5 


lOllC 


0.399948-0.399966 


8 


9 


lOllOlOlC 


0.402190-0.402194 


9 


7 


lOllOlC 


0.403530-0.403564 


01 


9 


lOllOlllC 


0.404690-0.404691 


00 


8 


lOllOllC 


0.406054 


12 


3 


IOC 


0.40912-0.41091 




2x3 


lOOlOC 


0.41092-0.41175 


13 


8 


lOOlOllC 


0.414432 



The forced Brusselator had been extensively studied with symbolic dynamics of ID maps*^. 
An incomplete U-sequence up to period six along the axis A — 0.46 — 0.2a; had already been 
found by Hao et. al'' which is also listed in Table 1. Only the periodic window lOOOlC was 
missing. Our investigation on the Poincare map with symbolic dynamics shows that L = 2 
for the parameter range 0.8056 < a; < 0.8194 so that the U-sequence up to period 6 might 
be incomplete. Recently, J.X. Liu has confirmed that the missing period lOOOlC is pruned^^. 



V. CONCLUSION AND DISCUSSION 

In this paper the universal bifurcation property and its hmitations of the Henon map in 2D 
parameter space (a, b) is discussed with symbolic dynamics. Two topological quantities 5 
and L are defined to characterize this topological universality. In the universal-L region, as 
that of ID unimodal map, there is a perfect MSS-sequence up to period<L along any axis 
provided that the axis is never tangent with any isoperiodic lines, though the phase portraits 
of the Henon map exhibit very complicated 2D behaviour. Extending this idea to many other 
two- or higher-dimensional systems ensures that the symbolic dynamics of ID mappings is an 
effective technique to investigate the imiversality in these two- or higher- dimensional systems 
and then the parameter for definite periodic motion may be predicted^^. We have presented 
two examples of ordinary differential equations (ODE's), the Rossler equations^^ and the 
forced Brusselator^, to demonstrate the validity and robustness of our approach. 
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It should be noted that only the short period is considered although the theory presented 
in this paper is also valid for higher period. In fact, in real experiment (or numerical study 
on ODE's or PDE's), only short periodic orbit can be obtained. Our investigation shows 
that it is not a surprizing result that universal sequences with short period are found in 
many experiments. Moreover, our result shows that it is also a universal phenomenon that 
universal sequences with short period can be found in many nonlinear dissipative systems. 
This observation ensures that the parameter of many periodic motion for many dynamical 
systems (such as some fluid system, e. g. ref. 23) can be well predicted. 

In this paper wc also show that even on a axis out of the universal region, the Henon map 
can exhibit approximately ID behaviour. This observation interprete the numerical results 
that in some nonlinear dynamical systems only incomplete U-sequences had been found^. 
Anyway, our defined universal-M region gives the background to interprete the experimental 
and numerical observations that complete or imcomplete U-sequences with short period can 
be found in many dissipative systems, and understand the limitations that ID symbolic 
dynamics can be used to study two- or high-dimensional dissipative systems. 
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Figure Captions 

Fig. 1 The Logistic map Xn+i = 1 — A*a;^ with jj, — 1.75488 exhibits a 3-cycle of the type 
IOC. 

Fig. 2 The strange attractor of the Henon map for {a,b)= (1.4, 0.16). The heavy lines 
outline the strange attractor. The diamonds are the "primary" tangencies. The dotted 
line connecting them divides the full attractor into two subsets marked by and 1. 

Fig. 3 The symbolic plane of the Henon map for {a,b)= (1.4, 0.16). (b) an enlarged part 
of the symbolic plane. 
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Fig. 4 The isoperiodic lines together with the universal-7 region (the region between 
two heavy hnes) and the borders for the Henon map exhibit an attracting set with 
initial points (a;o,yo) very close to original point (0,0) (the two heavy short dashed 
lines). The long dashed (1) and dash-dotted (11) lines represent two axes in and out 
of the universal-7 region, on which a complete and incomplete U-sequence is found 
respectively, (b) an enlarged part. 
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